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Abstract 

Weighted voting games (WVG) are coalitional games in which an agent's contribution to a 
coalition is given by his weight, and a coalition wins if its total weight meets or exceeds a 
given quota. These games model decision-making in political bodies as well as collaboration 
and surplus division in multiagent domains. The computational complexity of various solution 
concepts for weighted voting games received a lot of attention in recent years. In particular, 
Elkind et al.(2007) studied the complexity of stability-related solution concepts in WVGs, 
namely, of the core, the least core, and the nucleolus. While they have completely characterized 
the algorithmic complexity of the core and the least core, for the nucleolus they have only 
provided an NP-hardness result. In this paper, we solve an open problem posed by Elkind et 
al. by showing that the nucleolus of WVGs, and, more generally, fc-vector weighted voting 
games with fixed k, can be computed in pseudopolynomial time, i.e., there exists an algorithm 
that correctly computes the nucleolus and runs in time polynomial in the number of players n 
and the maximum weight W. In doing so, we propose a general framework for computing the 
nucleolus, which may be applicable to a wider of class of games. 

1 Introduction 

Both in human societies and in multi-agent systems, there are many situations where individual 
agents can achieve their goals more efficiently (or at all) by working together. This type of scenarios 
is studied by coalitional game theory, which provides tools to decide which teams of agents will 
form and how they will divide the resulting profit. In general, to describe a coalitional game, one 
has to specify the payoff available to every team, i.e., every possible subset of agents. The size 
of such representation is exponential in the number of agents, and therefore working with a game 
given in such form is computationally intensive. For this reason, a lot of research effort has been 
spent on identifying and studying classes of coalitional games that correspond to rich and practically 
interesting classes of problems and yet have a compact representation. 

One such class of coalitional games is weighted voting games, in which an agent's contribution 
to a coalition is given by his weight, and a coalition has value 1 if its total weight meets or exceeds a 
given quota, and otherwise. These games model decision-making in political bodies, where agents 
correspond to political parties and the weight of each party is the number of its supporters, as well 
as task allocation in multi-agent systems, where the weight of each agent is the amount of resources 
it brings to the table and the quota is the total amount of resources needed to execute a task. 

An important issue in coalitional games is surplus division, i.e., distributing the value of the 
resulting coalition between its members in a manner that encourages cooperation. In particular, it 
may be desirable that all agents work together, i.e., form the grand coalition. In this case, a natural 
goal is to distribute the payoff of the grand coalition so that it remains stable, i.e., so as to minimize 
the incentive for groups of agents to deviate and form coalitions of their own. Formally, this intuition 
is captured by several related solution concepts, such as the core, the least core, and the nucleolus. 
Without going into the technical details of their definitions (see Section[3]l, the nucleolus is, in some 
sense, the most stable payoff allocation scheme, and as such it is particularly desirable when the 
stability of the grand coalition is important. 

The stability-related solution concepts for WVGs have been studied from computational per- 
spective in Q. There, the authors show that while computing the core is easy, finding the least core 
and the nucleolus is NP-hard. These computational hardness results rely on all weights being given 
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in binary, which suggests that these problems may be easier for polynomially bounded weights. In- 
deed, paper [5] provides a pseudopolynomial time algorithm (i.e., an algorithm whose running time 
is polynomial in the number of players n and the maximal weight W) for the least core. However, 
an analogous question for the nucleolus has been left open. 

In this paper, we answer this question in affirmative by presenting a pseudopolynomial time 
algorithm for computing the nucleolus. 

Theorem 1. For a WVG specified by integer weights w±, . . . , w n and a quota q, there exists a 
procedure that computes its nucleolus in time polynomial in n and W = max.; Wi. 

As in many practical scenarios (such as e.g., decision-making in political bodies) the weights 
are likely to be not too large, this provides a viable algorithmic solution to the problem of find- 
ing the nucleolus. Our approach relies on solving successive exponential-sized linear programs by 
constructing dynamic-programming based separation oracles, a technique that may prove useful in 
other applications. 

A proof of Theorem[T]is presented in Section[4] after preliminaries in Section[3]and a discussion 
of related work in Section|2] The text rounds up by Section[5]that discusses conclusions and future 
work directions. 

2 Related work 

Another approach to payoff distribution in weighted voting games is based on fairness, i.e., dividing 
the payoff in a manner that is proportional to the agent's influence. The most popular solution 
concepts used in this context are the Shapley-Shubik power index |16| and the Banzhaf power 
index [ 1 1. Both of these indices are known to be computationally hard for large weights IT3l l4l. yet 
efficiently computable for polynomially bounded weights ifTUl . 

The concept of the nucleolus was introduced by Schmeidler lfl4l in 1969. Paper lfl4l explains 
how the nucleolus arises naturally as "the most stable" payoff division scheme, and proves that the 
nucleolus is well-defined for any coalitional game and is unique. Kopelowitz |9| proposes to com- 
pute the nucleolus by solving a sequence of linear programs; we use this approach in our algorithm. 

The computational complexity of the nucleolus has been studied for many classes of games, such 
as flow games 1 3 1 , cyclic permutation games [ 1 7 1 , assignment games 1 1 2 1 , matching games [ 8 1 , and 
neighbor games [7], as well as several others. While some of these papers provide polynomial-time 
algorithms for computing the nucleolus, others contain NP-hardness results. 

The work in this paper is inspired by [5 1, which shows that the least core and the nucleolus 
of a weighted voting game are NP-hard to compute. It also proves that the nucleolus cannot be 
approximated within any constant factor. On the positive side, it provides a pseudopolynomial time 
algorithm for computing the least core, i.e., an algorithm whose running time is polynomial in n 
and W (rather than in the game representation size 0(n log W)), as well as a fully polynomial time 
approximation scheme (FPTAS) for the least core. However, for the nucleolus, paper 13 contains 
no algorithmic results. 

3 Preliminaries and Notation 

A coalitional game G — (I, v) is given by a set of agents / = {1, . . . , n} and a function v : 2 1 — > K 
that maps any subset (coalition) of the agents to a real value. This value is the total utility these 
agents can guarantee to themselves when working together. A coalitional game is called simple if 
v(S) G {0, 1} for any coalition S C /. In a simple game, a coalition S is called winning if v(S) = 1 
and losing otherwise. 

A weighted voting game is a simple coalitional game G given by a set of agents I = {1, . . . , n}, 
their non-negative weights w = (wi ,...,«)„), and a quota q; we write G = (I; w; q). As the focus 
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of this paper is computational complexity of such games, it is important to specify how the game is 
represented. In what follows, we assume that the weights and the quota are integers given in binary. 
This does not restrict the class of WVGs that we can work with, as any weighted voting game has 
such a representation [11]. 

For a coalition SCI, its value v(S) is 1 (i.e., S is winning) if J2i<=s Wi — 9 ; otherwise, 
v(S) = 0. Without loss of generality, we assume that the value of the grand coalition / is 1, that is, 
Wi — 1' Also, we set W = max^j W{. 

For a coalitional game G = (I, v), an imputation is a vector of non-negative numbers p = 
(pi, . . . ,p n ), one for each agent in /, such that J^ieiPi = K-0- We refer to pi as the payoff of 
agent i. We write p(S) to denote Y^iesPi- Similarly, w(S) denotes J2ieS w i- 

An important notion in coalitional games is that of stability: intuitively, a payoff vector should 
distribute the gains of the grand coalition in such a way that no group of agents has an incentive to 
deviate and form a coalition of their own. This intuition is captured by the notion of the core: the 
core of a game G is the set of all imputations p such that 

p(S) > v(S) for all S C I. (1) 

While the core is an appealing solution concept, it is very demanding: indeed, for many games of 
interest, the core is empty. In particular, it is well known that in simple games the core is empty 
unless there exists a veto player, i.e., a player that is present in all winning coalitions. Clearly, this 
is not always the case in weighted voting games, and a weaker solution concept is needed. 

We can relax the notion of the core by allowing a small error in the inequalities ((TJ. This 
leads to the notion of e-core: the e-core of a game G is the set of all imputations p such that 
p(S) > ^(S) — e for all SCI. Under an imputation in the e-core, the deficit of any coalition S, 
i.e., the difference v(S) — p(S) between its value and the payoff that it gets, is at most e. Observe 
that if e is large enough, e.g., e > 1, then the e-core is guaranteed to be non-empty. Therefore, a 
natural goal is to identify the smallest value of e such that the e-core is non-empty, i.e., to minimize 
the error introduced by relaxing the inequalities in ((TJ. This is captured by the concept of the 
least core, defined as the smallest non-empty e-core of the game. More formally, consider the set 
{e | e < 1, e-core of G is non-empty}. It is easy to see that this set is compact, so it has a minimal 
element E\. The least core of G is its ei-core. The imputations in least core distribute the payoff in 
a way that minimizes the incentive to deviate: under any p in the least core, no coalition can gain 
more than E\ by deviating, and for any e' < E\, there is no way to distribute the payoffs so that 
the deficit of every coalition is at most e'. However, while the least core minimizes the worst-case 
deficit, it does not attempt to minimize the number of coalitions that experience the worst deficit, 
i.e., £i, nor does it try to minimize the second-worst deficit, etc. The nucleolus is a refinement of the 
least core that takes into account these higher-order effects. 

Recall that the deficit of a coalition S under an imputation p is given by d(p, S) — v(S) — p(S). 
The deficit vector of p is the vector d(p) = (d(p, Si), ... , d(p, S^™)). where Si, ... , S2™ is a list of 
all subsets of / ordered so that d(p, Si) > d(p, S2) > • ■ • > d(p, S2™). I n other words, the deficit 
vector lists the deficits of all coalitions from the largest to the smallest (which may be negative). The 
nucleolus is an imputation 77 = (171, ... , r/ n ) that satisfies d(f]) <\ cx d(x) for any other imputation 
x, where <i cx is the lexicographic order. It is known [ 14] that the nucleolus is well-defined (i.e., an 
imputation with a lexicographically minimal deficit vector always exists) and is unique. 

4 Algorithm 

The description of our algorithm is structured as follows. We use the idea of |9|, which explains 
how to compute the nucleolus by solving a sequence of (exponential-size) linear programs. In 
Section |4~T1 we present the approach of [9|, and argue that it correctly computes the nucleolus. 
This material is not new, and is presented here for completeness. In Section 14.21 we show how 
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to design separation oracles for the linear programs in this sequence so as to solve them by the 
ellipsoid method. While a naive implementation of these separation oracles would require storing 
exponentially many constraints, we show how to replace explicit enumeration of these constraints 
with a counting subroutine, while preserving the correctness of the algorithm. The arguments in 
Sections 14. II and l4~2l apply to any coalitional game rather than just weighted voting games. 

In Section l4~3l we show that for weighted voting games with polynomially-bounded weights the 
counting subroutine used by the algorithm of Section l4~2l can be efficiently implemented. Finally, 
in Section l4~4l we show how to modify this subroutine to efficiently identify a violated constraint if 
a given candidate solution is infeasible. The results in Sections [4. 3l and l4~4l are specific to weighted 
voting games with polynomially bounded weights. 

4.1 Computing the nucleolus by solving successive linear programs 

As argued in (5], the nucleolus can be computed by solving at most n successive linear programs. 
The first linear program CP contains the inequality p(S) > v(S) — e for each coalition SCI, 
and attempts to minimize e subject to these inequalities, i.e., it computes a payoff in the least core 
as well as the value e 1 of the least core. Given a (relative) interior optimizer (p^e 1 ) for CP 1 
(i.e., an optimal solution that minimizes the number of tight constraints), let E 1 be the set of all 
inequalities in CP 1 that have been made tight by p 1 (we will abuse notation and use E 1 to refer 
both to these inequalities and the corresponding coalitions). We construct the second linear program 
CP 2 by replacing all inequalities in E 1 with equations of the form p{S) = v(S) — e 1 , and try to 
minimize e subject to this new set of constraints. This results in e 2 < e 1 and a payoff vector p 2 
that satisfies p 2 (S') = v(S) - e 1 for all S G E 1 , p 2 (5) > v{S) - e 2 for all S E 1 . We repeat 
this process until the payoffs to all coalitions are determined, i.e., the solution space of the current 
linear program consists of a single point. It has been shown |9| that this will happen after at most n 
iterations: indeed, each iteration reduces the dimension of the solution space by at least 1. 

More formally, the sequence of linear programs {CP , . . . , CV n ) is defined as follows. The first 
linear program CP 1 is given by 



Let (p 1 , e 1 ) be an interior optimizer to this linear program. Let E 1 be the set of tight constraints 
for (p 1 , e 1 ) (and, by a slight abuse of notation, the coalitions that correspond to them), i.e., for any 
S e E 1 we have p 1 (S) = 1 - e 1 . 

Now, suppose that we have defined the first j — 1 linear programs CP 1 , . . . , CP 3 ~ l . For k = 
- 1, let (p k ,e k ) be an interior optimizer for CP k and let E fc = {S \ p k {S) = v(S) - e k }. 
Then the jth linear program CP 3 is given by 



min e subject to 




^Pi > v(S) -eforaU SCI. 



for all i = 1 , . . . , n 



(2) 




for all i = 1, . . . , n 



v(S) - e 1 for all SeE 1 



min e subject to < 



(3) 



Pi = v(S) - £ J ' _1 for all S C £* -1 



^p l >v{S)-e for all S U j k Z\ E fc . 
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Fix the minimal value of t such that there is no interior solution to CP 1 . It is not hard to see that 
the (unique) solution to LP 1 is indeed the nucleolus. Indeed, the nucleolus is a payoff vector that 
produces the lexicographically maximal deficit vector. This means that it: 

(i) minimizes e 1 such that all coalitions receive at least 1 — e 1 ; 

(ii) given (i), minimizes the number of coalitions that receive 1 — e 1 ; 

(iii) given (i) and (ii), minimizes e 2 such that all coalitions except for those receiving 1 — e 1 receive 
at least 1 — e 2 ; 

(iv) given (i), (ii) and (iii), minimizes the number of coalitions that receive 1 — e 2 , etc. 

Our sequence of linear programs finds a payoff vector that satisfies all these conditions; in particular, 
(ii) and (iv) (and analogous conditions at subsequent steps) are satisfied, since at each step we choose 
an interior optimizer for the corresponding linear program. The only issue that we have to address 
is that our procedure selects an arbitrary interior optimizer to the current linear program in order 
to construct the set E J . Conceivably, this may have an impact on the final solution: if two interior 
optimizers to CP 3 lead to two different sets TP , they may also result in different values of e 1+1 , 
so one would have to worry about choosing the right interior optimizer. Fortunately, this is not the 
case, as shown by the following lemma. 

Lemma 2. Any two interior optimizers (p, e) and (q, s) for the linear program CP 3 have the same 
set of tight constraints, i.e., the set T 3 is independent of the choice of the interior optimizer. 

Proof. First note that the set of all interior optimizers for CP is convex. Now, suppose that p and q 
are two interior optimizers for LP 3 , but have different sets of tight constraints. Then, by convexity, 
any convex combination ap + (1 — a)q of p and q is also an interior optimizer for CP 3 . However, 
the set of constraints that are tight for ap + (1 — a)q is the intersection of the corresponding sets 
for p and q, i.e., ap + (1 — a)q has strictly fewer tight constraints than p or q, a contradiction with 
p and q being interior optimizers for CP 3 . □ 

We conclude that when this algorithm terminates, the output is indeed the nucleolus. Next, we 
discuss how to solve each of the linear programs CP 3 , j = 1, . . . ,t. 

4.2 Solving the linear programs CV 1 , . . , LV 1 

It is well-known (see e.g. |fT31 |6j) that a linear program can be solved in polynomial time by the 
ellipsoid method as long as it has a polynomial-time separation oracle, i.e., an algorithm that, given 
a candidate feasible solution, either confirms that it is feasible or outputs a violated constraint. 
Moreover, the ellipsoid method can also be used to find an interior optimizer (rather than an arbitrary 
optimal solution) in polynomial time JU Thm. 6.5.5], as well as to decide whether one exists [6 
Thm. 6.5.6]. We will now construct a polynomial-time separation oracle for jth linear program LP 3 
in our sequence. 

It is easy to construct the part responsible for checking equations of LP 3 in (01, assuming that 
we already have an oracle for the (j — l)st program £P 3 ~ . Indeed, the latter oracle can be easily 
modified to also check whether the equation e = e 3 ^ 1 holds, thus providing an oracle for the optimal 
face of CP 3 . Then by (6, Thm. 6.5.5] we can compute a basis of the optimal face (which consists 
of at most n equations) in polynomial time. The separation oracle can then reject a candidate solution 
(p, e) if p violates one of those basis equations. 

Dealing with the inequalities of CP 3 is more complicated. A naive separation oracle would 
have to explicitly list the sets E 1 , . . . , TP~ l , which may be superpolynomial in size. Alternatively, 
one can treat this part of the oracle as a 0-1 integer linear feasibility problem, with 0-1 variables 
Xi encoding a set S £ U 3 k ~l i 1 Y, k that provides a separating inequality for the oracle input (p, e). 
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Namely, suppose that we have verified that p satisfies all the equations in CP 3 (as described above). 
Then, given an interior optimizer (p 3-1 , e^~ l ) for CP 3 ~ l , the valuesxi, . . . ,x n can be obtained as 
a solution to the following inequalities: 

5>r^ > (4) 

i 

^2piXi<l-e, (5) 

i 

y^^jXj > q. (6) 

i 

The problem with this approach is that for arbitrary rational p and p J 1-1 this system of inequali- 
ties is at least as hard as KNAPSACK, which is NP-complete. Moreover, as p and p 3 ^ 1 are pro- 
duced by the ellipsoid method, there is no guarantee that their bitsizes are small enough to use a 
(pseudo)polynomial-time algorithm for KNAPSACK. The only hope is to replace at least one of (HJl 
and (0 by something "tame". 

We will now present a more sophisticated approach to identifying a violated constraint. In a way, 
it can be seen as replacing checking with counting. Our construction proceeds by induction: to 
construct a separation oracle for CP 3 , we assume that we have constructed an oracle for CV° , 
and are given the sizes of sets S 1 , . . . , as well as the sequence (e 1 , . . . , e 3 ^ 1 ). 

By construction, any optimal solution (p, e) to CP 3 satisfies e < e J ' — 1 , so we can add the 
constraint e < e 3 ^ 1 to CP 3 without changing the set of solutions. From now on, we will assume 
that CP 3 includes this constraint. Our separation oracle will first check whether a given candidate 
solution (p, e) satisfies e < e 3 ^ 1 , as well as constraints Pi > for alii = 1, . . . , n and p(I) = 1, 
and reject (p, e) and output a violated constraint if this is not the case. Therefore, in what follows 
we assume that (p, e) satisfies all these easy-to-identify constraints. 

Now, a candidate solution (p, e) is feasible for CP 3 if p{S) = v(S) — e* for S 6 £*, t = 
1, . . . ,j - 1, and p(S) > u(S) - e for all S g ufZ^E*. Recall that the deficit of a coalition S C I 
under a payoff vector p is given by v(S) — p(S). Suppose that we have a procedure P(p, e) that, 
given a candidate solution (p, e), can efficiently compute the top j distinct deficits under p, i.e., 

m 1 = max{d(S) | S C 1} 

m 2 = m&x{d{S) | 5 C I, d(S) ^ m 1 } 

m? = max{d(S') | S C I, d{S) / m 1 , . . . , m^ 1 } 

as well as the numbers n 1 , . . . , n 3 of coalitions that have deficits of m 1 , . . . , m 3 , respectively: 

n k = \{S\S<^I 1 d(S) = m k }\, k = l,...,j. 

Suppose also that we are given the values e 1 , . . . , e 3 ^ 1 and the sizes s* of the sets £*, t = 1, . . . , j — 
1. 

Now, our algorithm works as follows. Given a candidate solution (p, e), it runs P(p, e) to obtain 
m t ,n t ,t= 1, . . . ,j. If e < e- 7-1 , it then checks whether 

(a) m* = e* and n l = s* for all t = 1 , . . . , j — 1 

(b) m 3 < e. 

If e = e- ?_1 , it simply checks whether m* = e* for t — I, ... ,j — I and n l = s* for all t = 
1, . . . , j — 20. If these conditions are satisfied, the algorithm answers that (p, e) is indeed a feasible 

1 Alternatively, if e = e J , one can verify whether (p, e) is a feasible solution to the previous linear program CP 3 ~ 1 . 
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solution, and otherwise it identifies and outputs a violated constraint (for details of this step, see 
Section |4~4l i. We will now show that this algorithm implements a separation oracle for CP 3 correctly 
and efficiently. 

Theorem 3. Given the values s t ,s t ,t— 1, ... ,j — 1, and a procedure V(p, e) that computes 
m t ,n t , t = 1, . . . , j, in polynomial time, our algorithm correctly decides whether a given pair 
(p, e) is feasible for CP 3 and runs in polynomial time. 

Proof. We start by proving an auxiliary lemma. 

Lemma 4. For any vector p and any t < j — 1 such that m s — e s , n s — s s for all s < t, the 
coalitions with deficit s s under p are exactly the ones in E s . 

Proof. The proof is by induction on s. For s = 1, we have that the largest deficit of any coalition 
under p is e 1 , and there are exactly s 1 coalitions with this deficit. Hence, (p, e 1 ) is an interior 
optimizer for CP , and therefore the lemma follows by Lemma|2] Now, suppose that the lemma has 
been proven for s — 1. By the induction hypothesis, p satisfies all constraints in E 1 , . . . , Also, 
under p there are at most s 1 + • • + s s_1 coalitions whose deficit exceeds e s , so for all coalitions 
not in U^z\T, r their deficit is at most e s . Finally, there are exactly s s coalitions whose deficit is 
exactly e s . Hence, (p,e s ) is an interior optimizer for LV S ', and therefore the lemma follows by 
Lemma |2] □ 

To prove the theorem, let us first consider the case e < e J_1 . Suppose that (p, e) satisfies (a) 
and (b). By using Lemma|4]with t = j — 1, we conclude that (p, e) satisfies all equations in CP 3 . 
Now, under p, m? is the largest deficit of a coalition not in Uj^E'. If this deficit is at most e, then 
the pair (p, e) is a feasible solution to CP 3 . 

Conversely, suppose that (a) or (b) is violated. If (p, e) satisfies (a) but violates (b), by us- 
ing Lemma [4] with t = j — 1 we conclude that, under p the coalitions in E* have deficit e* for 
t = 1, . . . ,j — 1, but the deficit of some coalition not in Uj~ 1 S* exceeds e. Hence, this coalition 
corresponds to a violated constraint. Now, suppose that (a) does not hold, and let s be the smallest 
index for which m s ^ e s or n s ^ s s . By using Lemma [4] with t = s — 1, we conclude that for 
r = 1, . . . , s — 1 the coalitions in E r have deficit e r . However, either the sth distinct deficit under 
p is not e s , in which case p violates a constraint in 2 1 \ \J s r l\Y7', or under p there are more than s s 
coalitions with deficit e s (note that by construction it cannot be the case that n s < s s ). In the latter 
case, there is a coalition in 2 1 \ U^ =1 E r whose deficit exceeds e s , thus violating the corresponding 
constraint. The case e = e 3 ^ 1 is similar. In this case for a candidate solution (p, e) to be feasible, 
it is not required that there are exactly s 3 coalitions with deficit e 3 ' -1 . Hence, the algorithm only 
has to decide if for alH = 1, . . . , j — 2, the coalitions with deficit e* under p are exactly the ones in 
E', and all other coalitions get at least e 3 ' 1 . Showing that our algorithm checks this correctly can 
be done similarly to the previous case. 

The bound on the running time is obvious from the description of the algorithm. □ 

To provide the value s 3 = |E J | for the subsequent linear programs CP 3+1 , . . . , CP n , we need 
to find an interior optimizer for CP 3 . Thm. 6.5.5 in (5) explains how to do this given a separation 
oracle for the optimal face, i.e., the set of all optimizers of CP 3 . Observe that such an oracle can 
be obtained by a slight modification of the oracle described above. Indeed, the optimal face is the 
set of solutions to the linear feasibility problem given by the constraints in CP 3 together with the 
constraint e — e 3 . The modified oracle first checks the latter constraint, reports the violation (and the 
corresponding inequality) if it happens, and otherwise continues as the original oracle. Clearly, the 
modified oracle runs in polynomial time whenever the original one does. Hence, we can compute 
s 3 in polynomial time by computing an interior solution (p, e) to CP 3 according to [6 Thm. 6.5.5], 
running P(p, e) to find n 3 , and setting s 3 = n 3 . 



7 



4.3 Implementing the counting 

We will now show how to implement the counting procedure V(p, e) used in Section |4~2l for WVGs. 
The running time of our procedure is polynomial in the number of players n and the maximum 
weight W. 

Our approach is based on dynamic programming. Fix a WVG (I; w; q), a payoff vector p, and 
j < n. For all k = 1, . . . , n, w = 1, . . . , nW, let X k w , ■■ ■ , X k be the bottom j distinct payoffs 
to coalitions in {1, . . . , k} of weight w, i.e., define 

x k,w = min{p(S') | S C {1, . . . , k}, w(S) = w} 

X\ w = min{p(S) | S C {1,. .., k},w(S) = w,p(S) + X\ >w \ 

X{ w = min{p(5) | S C {1, ... , k},w(S) = w,p(S) + X\^ . . . , X^} 
and let Y kw , . . . , Y kw be the numbers of coalitions that get these payoffs, i.e. set 

Yl w = \{S | S C {1,. ..,fc},ti;(5) = w,p(5) = fori = 1, . . . , j. 

These quantities can be computed inductively for k = 1, . . . , n as follows. 

For k — 1, we have X| ^ = pi if iy = w\ and +oo otherwise, ij 1 = 1 if w = w% and 
otherwise, and X\ w = +oo, Y{ w = for t = 2, . . . ,j. 

Now, suppose that we have computed x l_ ljW , x l-i yW ,Y^_ lw , . . . , Y^_ lw for all w = 
1, . . . , nW. Consider S C {1, . . . , fc} receiving one of the bottom j distinct payoffs to subsets of 
{!,..., k] of weight to, i.e., p{S) G {X^, ■ ■ • , Then either 

(1) "S 1 C {1, . . . , k — 1}, in which case 5 must be among the coalitions that receive one of the 
bottom j distinct payoffs to subsets of {1, . . . , k — 1} of weight w, i.e., we have p(S) £ 

(2) k G S, in which case S \ {k} must be among the coalitions that receive one of the bottom j 
distinct payoffs to subsets of {1, .... k — 1} of weight w — Wk, i.e., we have p(S \ {k}) G 

Consider the multi-set S k , w = { x l_ 1<w , ■ • .,X J k _ lw ,p k + Xj._ hw _ Wk ,. .., Pk + Xl_ lw _ w J. 
By the argument above, we have 

x k,w = min{a; | x G 

= min i x I a: G S k<w ,x=f Xl w } 

x l,w = m[n i x I x G s k,w,x=£ Xl wl ...,X 3 -^}. 

The number of coalitions that receive the payoff X l k w , i.e., Y kw , t = 1, . . . , j, depends on how 
many times X k appears in S kjW . If it only appears once, then there is only one source of sets 
that receive a payoff of X\ w , i.e., we set Y kw — Y k s _ 1 w if X\ w appears as X k-1 w for some 
s = 1, . . . , j, and we set Y^ w = Y k s _ l w _ Wk if X\ w appears as X%_ hw _ Wk + p k for some s = 
1, . . . , j. On the other hand, if X k w appears twice in S kyW (first time as X k _ 1 w and second time as 
p k + X k _ 1 w _ Wk for some s,s' — 1, . . . , j), we have to add up the corresponding counts, i.e., we 
set V* — Y s A- Y s ' 

After all X\ w , ■ ■ ■ , X^ w , Y^ w , . . . , Y£ ' w have been evaluated, it is not hard to compute m f , n f , 
t = 1, . . . , j. Indeed, the top j deficits appear in the multi-set 

S = {l w - X\ w , ...,I W - X 3 nw I to = 1, ... , nW}, 
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where I w — 1 if w > q and I w = if w < q (recall that q is the quota of the game, i.e., v(S) = 1 if 
and only if w(S) > q). Hence, we can set 

m = max{i | x G S} 

m = max{i | x £ S, x ^ m 1 } 

m? = max{s | x s <S, x ^ m , . . . , m-* -1 }. 

The procedure for computing n 1 , t — 1, . . . , j, is similar to that of computing Y]£ (see above): we 
have to check how many times m* appears in S and add the corresponding counts. 

In the next subsection, we will show how to find a violated inequality if (p, e) is not a feasible 
solution to CP 3 . 

4.4 Identifying a violated constraint 

Consider CP 3 and a candidate solution (p, e). Suppose that the algorithm described in the previous 
subsection has decided that (p, e) is not a feasible solution to CP 3 . This can happen in three possible 
ways. 



(a) 


m s 


= e s , 


n s 


= s s for s = 


1,.. 


..,£- 


l,but 


m e 7^ e e for an £ < j. 


(b) 


m s 


= e s , 


n s 


= s s for s = 


1,. 


..,£- 


1, m 


— e , but n e ^ s e for an £ < j 


(c) 


m s 


= e s 


n s 


= s s for s = 


1,.. 


■J- 


l,but 


m? > e. 



In cases (a) and (b), there is a violated equation in (0, while in (c) there are none (but there is a 
violated inequality). Thus (a) and (b) can be handled using the ideas discussed in the beginning of 
Section l4!2l Indeed, as argued there, we can efficiently compute the basis of the optimal face of the 
feasible set of CP 3 ~ using the ellipsoid method. One can then easily check if a candidate solution 
violates one of the equations in the basis (recall that there are at most n of them), and, if this is the 
case, report one that is violated. Hence, we only need to show how to identify a violated constraint 
in case (c). However, for completeness, we present here a purely counting-based algorithm for each 
of the cases. 

In case (a), let (p, e l ) be an interior optimizer for CP 1 . Under p, the deficit of any coalition 
in 2 1 \ uijZ^X" is at most e e . On the other hand, under p, there are n l coalitions in 2 1 \ U^Z^S 3 
whose deficit is m 1 > e . Each of these coalitions corresponds to a violated constraint: indeed, if 
such a coalition is in S s , s > £, then CP 3 requires that its deficit is e s < e e < m e , and if it is in 
2 / \U^£ s ,then CP 3 requires that its deficit is at most e < e < m e , Hence, it suffices to identify 
a coalition whose deficit under p is m l . To this end, we can modify the dynamic program for p as 
follows. Together with every variable Xj, w , t = 1, . . . , j, k — 1, . . . , n, w = 1,..., nW, we will 
use an auxiliary variable Zt which stores a coalition whose payoff under p is equal to X\. w . The 
values of Z\ can be easily computed by induction: if X\ = X^._ 1 for some s = 1, . . . , j then 
= Z k-i, w >^XXk, w = Pk + XU, w - Wk for some a = 1, . . . then Z{ w =^|_ 1)W U{fc} 
( if Xl,w = X k-i, w = Pk + X k-i, w ~w k ' we can set Z k.w to e i tner of these values). Now, there exist 
some t and w such that w > q and 1 — X^ w = mr or w < q and — w = m , such t and w can 
be found by scanning all X l n w . The corresponding set Z l n w has deficit m e under p, and therefore 
corresponds to a violated constraint. 

In case (b), as before, let (p, e) be an interior optimizer to CP 3 ~ l . There exists a coalition whose 
deficit under p is e , but whose deficit under p is strictly less than e . To find such a coalition, run 
P(p, e) in order to compute the corresponding values X\, , Y£ t = 1, . . . , j — 1, k = 1, . . . , n, 
w = 1, . . . , nW. Define Z w as follows: if there exists some t G {1, . . . , j — 1} such that X^ w = 
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I w — e l , set Z w — Y* ', otherwise, set Z w = 0. Z w is defined similarly: if there exists an 
s S {1, . . . , j — 1} such that w = I w — er, set Z w = Y% w ; otherwise, set Z w = 0. The variables 
Z w and Z w count the number of coalitions that have total weight w and have deficit e under p 
and p, respectively. We have n e = J2w=i...., n w Z w ^ = Y, w =i,...,nW As n e > s e , there 
exists a weight w such that Z w > Z w . Set q = I w — e , i.e., q is the total payment received by the 
coalitions counted by Z w and Z w . Now, we have Z w — Z™ + Z~ n , where Z™ is the number of 
coalitions of weight w that include n, have weight w and receive total payment q, and Z~ n is the 
number of coalitions of weight w that do not include n, have weight w and receive total payment 
e , Z™ and Z~ n can be defined similarly. We can easily compute these quantities: for example, 
Z™ is the number of subsets of {1, . . . , n — 1} that have weight w — w n and receive total payment 
q - p n , i.e. Z™ = Y*_ 1>w _ Wn if there exists a t G {1, . . . , j - 1} such that X l n _ lw ^ Wn =q- p n , 
and Z~ n = otherwise. It follows immediately that Z™ > Z™ or Z^" > Z^" (or both), and we 
can easily verify which of these cases holds. In the former case, we can conclude that the number 
of coalitions in {1, . . . , n — 1} that have weight w — w n and are paid q — p n under p exceeds the 
number of coalitions in {1, . . . , n — 1} that have weight w — w n and are paid q — p n under p. In 
the latter case, we can conclude that the number of coalitions in {1, . . . , n — 1} that have weight w 
and are paid q under p exceeds the number of coalitions in {1, . . . , n — 1} that have weight w and 
are paid q under p. Continuing in the same manner for n — 1, .... 1, we can identify a coalition that 
is paid q under p, but not under p. 

Case (c), i.e. m? > e, is similar to (a) and can be handled in the same manner. 

5 Conclusions and future work 

In this paper, we proposed a new technique for computing the nucleolus of coalitional games. 
Namely, we have shown that, when constructing the separation oracle for the jth linear program 
CP 3 , instead of storing the sets of tight constraints for the linear programs CP 1 , t = 1, . . . , j — 1, 
it suffices to store the sizes of these sets as well as the top j — 1 deficits of an interior optimizer 
(p J _1 , e) to CV J ~ . A feasibility of a candidate solution (p. e) to CP 3 can then be verified, roughly, 
by computing the top j deficits for p as well as the number of coalitions that have these deficits, and 
comparing these values to their pre-computed counterparts for (p J_1 , e). 

We then demonstrated the usefulness of this technique by showing that for weighted voting 
games with polynomially-bounded weights both the top j deficits and the number of coalitions that 
have these deficits can be efficiently computed using dynamic programming. This allows us to 
implement the separation oracles for our linear programs in pseudopolynomial time. Combining 
this with the ellipsoid algorithm results in a pseudopolynomial time algorithm for the nucleolus 
of weighted voting games, thus solving an open problem posed by [5|. Furthermore, the general 
technique put forward in this paper effectively reduces the computation of the nucleolus to solving 
a natural combinatorial problem for the underlying game. Namely, we can state the following meta- 
theorem: 

Theorem 5. Given a coalitional game G, suppose that we can, for any payoff vector p, identify 
the top n distinct deficits under p as well as the number of coalitions that have these deficits in 
polynomial time. Then we can compute the nucleolus of G in polynomial time. 

We believe that this framework can be useful for computing the nucleolus in other classes of 
games. Indeed, by stripping away most of the game-theoretic terminology, we may be able to find 
the nucleolus by applying existing results in combinatorics and discrete mathematics in a black-box 
fashion. 

In the context of weighted voting games, our assumption that the weights are polynomially 
bounded (or, equivalently, given in unary) is essential, as @ shows that the nucleolus is NP-hard to 
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compute for WVGs with weights given in binary. Moreover, in many practical scenarios the agents' 
weights cannot be too large (e.g., polynomial functions of n), in which case the running time of our 
algorithm is polynomial. 

By a slight modification of our algorithm, we can obtain a pseudopolynomial time algorithm for 
computing the nucleolus in k-vector weighted voting games for constant k. Informally speaking, 
these are games given by the intersection of k weighted voting games, i.e., a coalition is considered 
to be winning if it wins in each of the underlying games. There are some interesting games that can 
be represented as fc-vector weighted voting games for small values of k (i.e., k — 2 or k = 3), but not 
as weighted voting games, most notably, voting in the European Union jZ|. Hence, this extension of 
our algorithm enables us to compute the nucleolus in some real-life scenarios. The overall structure 
of our algorithm remains the same. However, the dynamic program has to be modified to keep track 
of several weight systems simultaneously. 

Another natural way to address the problem of computing the nucleolus is by focusing on ap- 
proximate solutions. Indeed, proposes a fully polynomial time approximation scheme (FPTAS) 
for several least-core related problems. It would be natural to expect a similar result to hold for 
the nucleolus. Unfortunately, this approach is ruled out by [5 1, which shows that it is NP-hard to 
decide whether the nucleolus payoff of any particular player is 0, and therefore approximating the 
nucleolus payoffs up to any constant factor is NP-hard. Nevertheless, one can attempt to find an 
additive approximation to the nucleolus, i.e., for a given error bound S > 0, find a vector x such that 
\f]i — Xi\ < S for i = 1, ... ,n. This can be useful in situations when the agents' weights cannot 
be assumed to be polynomially bounded with respect to n, e.g., in the multiagent settings where 
the weights correspond to agents' resources. We are currently investigating several approaches to 
designing additive approximation algorithms for the nucleolus. 
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